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Abstract. We present a generalization of Kruzkov's theory to manifolds. Nonlin- 
ear hyperbolic conservation laws are posed on a differential (n + l)-manifold, called a 
spacetime, and the flux fleld is deflned as a fleld of n-forms depending on a parameter. 
The entropy inequalities take a particularly simple form as the exterior derivative of a 
family of n-form fields. Under a global hyperbolicity condition on the spacetime, which 
allows arbitrary topology for the spacelike hypcrsurfaccs of the foliation, we establish 
the existence and uniqueness of an entropy solution to the initial value problem, and we 
derive a geometric version of the standard semi-group property. We also discuss an 
alternative framework in which the flux field consists of a parametrized family of vector 
fields. 
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1. Introduction. Hyperbolic conservation laws on manifolds arise in many 
applications to continuum physics. Specifically, the shallow water equations for perfect 
fluids posed on the sphere with prescribed topography provide an important model 
describing geophysical flows and exhibiting complex wave structures (Rosby waves, etc.) 
with several distinct scales. In [J— [5] and 1131 1181 1191 120| , the author together with 
several collaborators recently initiated a research program and posed the foundations for 
the analysis and numerical approximation of hyperbolic conservation laws on manifolds. 
These equations form a simplified, yet challenging, mathematical model for studying 
the propagation of nonlinear waves, including shock waves, and their interplay with the 
background geometry. 

2. Hyperbolic conservation laws based on vector fields. The 

standard shock wave theory posed on the Euclidian space covers nonlinear hyperbolic 
equations of the form 

n 

dtu + ^djf^(u) =0, u = u{t,x) G R, t > 0, X e K", 

where : R — > K are given functions, and, more generally, applies to balance laws with 
variable coefficients and source-terms. Pioneering works by Lax, Oleinik, and others 
in the 50's and 60's set the foundations about discontinuous entropy solutions to such 
equations. Later, well-posedness theorems within the class of entropy solutions were 
established by Conway and SmoUer (1966), Volpert (1967) (for solutions with bounded 
variation), Kruzkov (1970) (for bounded solutions), and finally DiPerna (1984) (for 
measure- valued solutions). Based on this mathematical theory, active research followed 
in the 80's and 90's that led to the development of robust and high-order accurate, 
shock-capturing schemes. 
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The generalization of the above theory to manifolds has only recently received par- 
ticular attention, driven by problems arising in geophysical science and general relativity. 
In the present section, we consider a class of conservation laws based on parametrized 
vector fields, and we suppose that M is a compact, oriented n-dimensional manifold, 
endowed with an L°° volume form such that, in a finite atlas of local coordinate charts 
(x^) and for uniform constants c < C, 

uj = LJdx^ ■ ■ ■ dx"- , 0<c<oj<C. 

Recall that the divergence of an L°° vector field X is defined in the sense of distributions 

by 

(div^ X,e) := f {de){X)u 
J M 

for all test-function 6» : Af -s> R. 

Definition 2.1. To any parameterized family of (smooth) vector fields f = fp{u) G 
TpM depending upon il G M, one associates the hyperbolic conservation law on the 
manifold (M, ui) 

(2.1) dtu + div^ (/(«)) =0, u = u(t, p), t>0, p€ M. 

The parametrized flux field f is called geometry-compatible if constants are trivial 
solutions, that is 

(2.2) {div^ f)p{u) = 0, ueR,peM 

in the distribution sense, A convex entropy pair is a convex function f/ : R — > R 
together with a parametrized family of vector field 

duU{v)dufp{v)dv GTpM, « G R, p e M. 

The equation I I2.1I I is geometric in nature and does not depend on a particular choice 
of local coordinates on M. The flux field is a section of the tangent bundle TM and, 
in general, there is no concept of "spatially constant flux" , unlike in the Euclidian case 
where one may arbitrarily choose a vector fp{u) at some point p G R" and parallel- 
transport it in R", trivially generating a constant (and, therefore, smooth) vector field. 
The geometry-compatibility condition is natural if one observes that it is satisfied by 
the shallow water equation. 

Remark 2.1. An important class of geometry-compatible conservation laws is based 
on projected gradient fiux fields, as it is called in J^, and is defined as follows. 
Consider the unit sphere embedded in the Euclidian space R^ and a function h : 
R'^ X R — R being fixed, and define 

fp{u) :=pA (VR3fe)p(?I) e TpS^, p G c R^ 

Now, to formulate the initial value problem associated with I I2.1I I we are given an 
initial data mq • Af — > R and we impose 

(2.3) m(0,p) = no(p), pGM. 

It is well-known that entropy inequalities must be imposed on weak solutions to con- 
servation laws. In this juncture, it is interesting to point out the following property of 
geometry-compatible flux field: If -u = u{t,p) is a smooth solution and U = U{u) is an 
arbitrary function, then clearly 

dtU(u) + duU{u) div^ {fp{u)) = 0, 

and it can then be checked that a necessary and sufficient condition for the existence of 
a vector field F = Fp{u) G TpM such that 

duU{v)div^ {fp{v)) = div„ {Fp{v)) 
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for all functions v = v{p) is precisely that the flux fp{u) be geometry-compatible. This 
observation motivates the following definition which impose the entropy inequality 

(2.4) dtU(u) + div^ (Fiu)) <0 

in the sense of distributions for all convex U. 

Definition 2.2. An entropy solution to the conservation law i2. IV . posed 
on the manifold (M, cj) and assuming the initial data uq g L°°(M), is a scalar field 
u G L°° ^R"*" X satisfying the entropy inequalities 

f I (u(u)dte + de{F{u)))ujdt+ f u(uo)e{o,-)i^>o 

Jr+ Jm ^ ' J M 

for all convex {U, F) and test-functions 6 > 0. 

Note that integration is performed w.r.t. the volume form on M, and the differential 
d6 is a 1-form field acting on vector fields. Our main existence result is as follows. 

Theorem 2.1 (Well-posedness theory on a manifold). Let {M,lj) be a compact 
n-manifold endowed with an L°° volume form, and f = fp{u) S TpAI be a geometry- 
compatible flux field. 

1. Given any data uq € L°°{M), the initial value problem associated with the 
conservation law \2.1\l admits a unique entropy solution u € L°°(R+ X M). 

2. Entropy solutions satisfy the following stability property for all q G [l,oo] 
and t2 > ti 

(2-5) ll"(*2)|li,9 (M) < (M)- 

3. Moreover, the following semi-group property holds for any two entropy 
solutions u,v and t2 > ti 

(2.6) \\v{t2) - M(t2)|Ll^(M) < - «(*l)ILi(A/)- 

This theorem provides a generalization of Kruzkov's theory to a manifold, endowed 
with a volume form with possibly low regularity. The stability and contraction 
properties provide geometry- independent bounds, which are useful to design approxi- 
mation schemes. For further discussions on the well-posedness theory based on vector 
fields, we refer the reader to Ben-Artzi and LeFloch 4 and Amorim, Bcn-Artzi and 
LeFloch [J, as well as Panov [211 [22] . 

Consider a geometry-compatible conservation law with flux field / together with 
some initial data no 6 L°°{M). Then, an entropy measure- valued solution to 
ll2ll-|[23t. by definition, is a measure-valued field {t,p) G X M »— > I't^p G Prob(R) 
such that 

/ / ({iyt,p,u)dte-\-{{ut.p,F),de)]w{p)dt-\- f u(uo)e(o,-)uj>o 

Jr+ Jm V / Jm 

for all convex entropy pairs and test-functions 8 = 0{t,p) > 0. Following DiPerna 
we can establish the following uniqueness result. 

Theorem 2.2 (Uniqueness and compactness framework). Ifu is an entropy measure- 
valued solution to the conservation law 12. l[l posed on the manifold (M, u)) and assuming 
some initial data uq £ L°°(M) at t = 0, then 

i^t.p = Su{t,p) a.e. p e M. 

where u S L°°(R"'" X M) denotes the unique entropy solution to the same problem. 

This theorem is very useful to establish the convergence of approximation schemes 
since, in comparison with the standard technique based on a total variation bound, it 
requires weaker and geometrically more natural bounds on (the entropy dissipation of) 
the approximate solutions, as well as weaJcer assumptions on the manifold geometry. 

In the Euclidian case, the total variation diminishing (TVD) property was found to 
be very useful for the development of high-order schemes. On manifolds, provided ui is 
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sufficiently smootii and tiic initial data uq has bounded variation then u{t) has bounded 
variation for all < > 0. In fact, tiiere exists C > depending on the geometry of M and 
II"oIIl°°(M) such that 

TV(u(t)) < e^^ {l+TV(uo)), 

and the total variation might be unbounded as t — > oo. 

Introduce the total variation along a given vector field X by 

(2.7) TVx{u):= sup / Mdiv„(0X)a;. 

Theorem 2.3 (Total variation estimate). Consider a geometry- compatible conser- 
vation law i2.1\) posed on the manifold {Al, lu) and suppose that the initial data have 
bounded variation. Then, the total variation diminishing property (TVD) 

TVxl,u(t2)) < TVxHti)), t2 > ti, 

holds provided the flux field f and the given field X satisfy the compatibility condition 

(2.8) Cx{duf{u)) ■.= [aufiu),x]=0. 

On the other hand, from the time-invariance of the equation and the contraction 
property, it follows that 

(2.9) WdiVu^UHmiM^iM) ■■= sup / de{f(u(t)))uj 

\\e\\^o<iJM 

is diminishing in time. 

The condition 12.81 1 is usually not satisfied in the examples of interest, except 
when there exists some decoupling into a one-parameter family of independent, one- 
dimensional equations. (See below.) Hence, the curved geometry restricts the class of 
flux enjoying the TVD property, and we can not expect to rely on a TVD property 
to develop approximation schemes. In contrast, in the Euclidian case no restriction is 
implied on (p-independent) flux and the TVD property always holds. 

Remark 2.2. Returning to the projected gradient flux field on C introduced 
in Remark \2.1\ and choosing a vector field X of the form 

Xp:=ph{SJ^ik)p, peM^, 

then the condition Cx (duf{u)^ = on holds if and only if the vectors dufp{u) 

and Xp are parallel, and the proportionality factor C{u,p) is constant along the integral 
curves of X : 

dufp{u) = dnC{u,p)Xp, X{C{u,p)) = 0, MgR, peM. 



3. Hyperbolic conservation laws based on diflferential forms. 

In this section based on LeFloch and Okutmustur 1191 120| . we consider hyperbolic con- 
servation laws posed on an {n -\- l)-dimensional differentiable manifold M, referred to 
as a spacetime. We emphasize that, in the framework now developed, no geometric 
structure is a priori imposed on M. 

Definition 3.1. The hyperbolic conservation law on the differentiable {n-\- 
l)-manifold M associated with a parametrized family of L°° n-form fields p a M t-^ 
i^p{u) depending smoothly upon the parameter u, by definition, reads 

(3.1) d(uj{u))=Q, 

whose unknown is the scalar field -u : M — > R. The parametrized flux field ui is called 
geometry-compatible if it is closed, that is, for each « G M, {du]){u) = in the 
distribution sense. 
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By definition, weak solutions u G L°°(M) must satisfy the following weak form of 
the conservation law 



(3.2) / d9 A uj{u) = 



for all test-functions 9 : M 

M 



Then, to pose the initial value problem we impose the following global hyperbolicity 
condition: M admits a smooth foliation by compact, oriented hypersurfaces, i.e.. 



M=\Jn 



t>0 



Furthermore, denoting by i-^^ : Ht — > M the injection map, we impose that each 
hypersurface is spacelike in the sense that the n-forms 

(3.3) dn(iH^io{u)^ > 

are (positive) n-volume forms for each t > and all relevant u. 

To formulate the initial value problem associated with 113.11 1 we are given some initial 
data uq G L°°(Wo) and we search for a scalar field u g L°°(M) satisfying 

(3.4) u^-Hg = UQ. 

In the present framework, a convex entropy fiux is now a parametrized family of 
n-forms Q = Clp{u) such that there exists a convex [7 : M — > M 



/u 
duU dui^p du' , p 



e M, u e : 



If the flux uj is geometry-compatible, then smooth solutions satisfy the additional con- 
servation laws 

d{n{u)) = 0, 

so that the entropy inequalities take the following form 

(3.5) d{n{u)) < 0. 

Definition 3.2. An entropy solution to the geometry-compatible conser- 
vation law is. 1\ } is a scalar field u = u{p) g L°°{M) satisfying for all convex entropy 
flux Q = Qp{u) and smooth functions 9 > 



f d9An{u)+f (j*n)(Mo) > 0. 

JM J-Hn 



Introduce also Kruzkov's entropy fiux fields, defined here as parametrized fields of 
n-forms: 



Q{u, v) := sgn (u — v) ^tj('(t) — '^(f)^ . 



Theorem 3.1 (Well-posedness theory on a spacetime. Geometry-compatible fiux). 
Consider a conservation law with a geometry-compatible flux u) = u){u), posed on a 
globally hyperbolic, spatially compact manifold M = Ut>0^*' Then, the initial value 
problem with data uq G L°°{'Ho) admits a unique entropy solution u G L°°{M), and 
for any convex entropy flux f2 the function t ^ J!^^ i'^^Q{u) is well-defined and non- 
increasing. Moreover, for any two entropy solutions u, v 

JHt 

is non-increasing in time. 

Theorem 3.2 (Well-posedness theory on a spacetime. General fiux). Consider 
a conservation law with general flux field u) = uj{u), posed on globally hyperbolic. 
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spatially compact M with (past) boundary Hq. Consider initial data uq satisfying 
fllg l*'Ho'^("o)| < oo- Then, there exists a semi-group of entropy solutions uq i—^ u = 
Suo and for all smooth functions 9 > 

(3.6) / dehQ.(u,v)+ iK„^^(MO,'^o) > 0. 
Jm J-Ho 

Moreover, for any two spacelike hypersurfaces "^2 in the future of Hi 

(3.7) / i^ fl{u,v) < / i^ fl{u,v) < oo. 
J-H2 JUi 



Remark 3.1. On the one- dimensional torus M = , the case of general flux 
fields was covered earlier by LeFloch and Nedelec \17^ under the assumption that f 
is strictly convex in u. Then, entropy solutions satisfy the generalized Lax formula 



u{t,x) = iduf) — v(*'^) ^ ^ which the function y = y{t,x) is determined by a 
minimization argument over the family of characteristics defined by 

a,X = (a„/o/-i)(^-^j, X(Q)=y, cG 



Interestingly, it was established in I17f that any two entropy solutions u, v : IR+ xT^ 
satisfy 

\\v{t2) - «(t2)|ll,l^{Tl) < lk(*l) - "(*l)lll,i,{Tl)> *2 > h. 



4. Intrinsic finite volume schemes. 

4.1. General intrinsic approach. Using the framework based on n-diffe- 
rential forms, Amorim, LcFloch, and Okutmustur [2] have introduced finite volume 
schemes, which generate piecewise constant approximations on unstructured triangu- 
lations. One considers here the initial value problem associated with the conservation 
law 113.11 1 and discretize the weak form 1 13.21 1 of the conservation law. Indeed, the fi- 
nite volume approach applies directly to the geometric weak form and, in particular, 
no local coordinates should be chosen at this stage. The finite volume discretization 
on a manifold M requires only the n-differential structure put forward in the previous 
section. 

One introduces a family of triangulations made of curved polyhedra K, each of them 
having a past boundary e^, a future boundary ej^, and "vertical boundary" elements 
e G dK^. For each boundary element e there is precisely one element i^^ such that 
K n Kg = e. Then, one defines the notion of "total flux" along spacelike hypersurfaces 

by setting 

/ i*a;(n) / i*ui(u ±) =: q ± (u ±), 

"k "k 

where n : M — > R denotes the entropy solution of the initial value problem under 
consideration. 

Then, applying Stokes theorem to the conservation law l|3.H l. we derive the finite 
volume scheme 

(4.1) g+(u+) = g (n )- V qe,x(«-,'"- ), 

eedKO 

in which the total flux along vertical faces are scalars determined from Lipschitz con- 
tinuous, numerical flux cie,K : — > R satisfying 
1. Consistency property: 

qX,e(«,M) = Qeiu), 
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2. Conservation property: 

3. Monotonicity property: 

9u<iK,e > 0, dv<lK,e < 

for all u,v £ M.. For instance, Lax-Friedrich flux or Godunov flux may be used to 
determine the functions qe,if • The use of total flux leads to a geometrically natural 
setting, in which the volume of elements K or boundary elements e do not arise, so that 
the proposed formulation is both conceptually and technically simpler. 

The convergence of finite volume schemes toward the entropy solution of the initial 
value problem is established in [21120], Let us mention here the main difficulties: possible 
blow-up of the TV norm; derivation of uniform estimates tied to the geometry, only 
involving n-dimensional volumes. The main issues in our analysis are as follows: 

1. Derivation of a discrete maximum principle. 

2. Derivation of discrete entropy inequalities. 

3. Derivation of entropy dissipation estimates. 

4. The sequence of finite volume approximation generate an entropy measure- 
valued solution Up which, by our generalization of DiPerna's uniqueness the- 
orem, coincides with i/p = 5„(p), which implies that the scheme converges 
strongly. 

DiPerna's measure-valued solutions were used to establish the convergence of sche- 
mes by Szepessy |23l 1241 . Coquel and LeFloch [8l l9l [T0|, and Cockburn, Coquel, and 
LeFloch |6l [7] . For many related results and a review about the convergence techniques 
for hyperbolic problems, we refer to Tadmor |25| and Tadmor, Rascle, and Bagnerini 
|26| . Further hyperbolic models, including also a coupling with elliptic equations, as 
well as many applications were successfully investigated by Kroners |14|, and Eymard, 
Gallouet, and Herbin |12| . For higher-order schemes, see the paper by Kroner, Noelle, 
and Rokyta 1151 . An alternative approach to the convergence of finite volume schemes 
was discovered by Westdickenberg and Noelle '27'. 

4.2. A second-order intrinsic scheme on the sphere. Ben-Artzi, 

Falcovitz, and LeFloch [5] adopted the fully intrinsic approach and extended the finite 
volume scheme to second-order accuracy (by using the generalized Riemann problem 
methodology), while implementing it for classes of hyperbolic conservation laws posed 
on the sphere S^. The sphere is the most important geometry for the applications to 
geophysical flows. Earlier works in this field were not fully geometric, and regarded the 
manifold as embedded in some Euclidian space and using projection techniques from 
the ambient space to the sphere. A fully intrinsic approach has definite advantages in 
terms of convergence and accuracy, and at the first-order at least, rigorous convergence 
analysis can be made. 

We present here some features of the proposed scheme and refer to [5] for further 
discussions, especially on the implementation. We consider the conservation law 

(4.2) dtu(t, x) + Vs2 ■ {F{x, u{t, x))) = 

with flux vector tangent to the sphere, in the general form 

F(x, u) = n{x) A ^{x, u). 

We are particularly interested in geometry-compatible flux vectors, and have proposed 
the following broad class of gradient flux vector 

^{x, u) = V/i(x, u). 

In particular, this includes homogeneous flux vectors having $ = $(«). Interestingly 
enough, even in this case, the corresponding flux covers a broad and non-trivial class of 
interest. 

Several classes of solutions of particular interest have been constructed, which exhibit 
very rich wave structure: 
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1. Spatially periodic solutions. 

2. Non-trivial steady state solutions. 

3. Confined solutions that remain compactly supported for all times. 
Explicit formulas for such solutions together with various choices of corresponding 

flux vector fields can be found in [5]. These formulas are useful for numerical investiga- 
tions. 

An intrinsic triangulation v^as used, a web-like mesh made of segments of longitude 
and latitude lines. Artificial coordinate-singularities at the poles are coped with by 
suitably adapting the mesh in a non-conformal fashion as one approaches the North 
and South poles. To derive explicit formulas for the flux, for instance, in the case of 
an homogeneous fiux vector 3>(?i), one may introduce the following basis of the tangent 
space 

ix = — sin A ii -f cos Xi2, = — sin 4> cos A ii — sin 4> sin A 42 + cos 4> is \ 

in terms of some Euclidian basis 11,12, i^irM? . Then, by setting 

= *l + /2('") *2 + h{u) 13, 

one finds 

F{x, u) = Fx(X, </>, u) ix + F^{X, </>, u) i^, 

with 

(4.3) Fx (A, </>, u) = fi (u) sin </> cos A -f /2(m) sin <f) sin A — f:j{u) cos ij>, 
and 

(4.4) F^{X,(p,u) = -fi{u) smX + f2{u) cosA. 

The intrinsic finite volume discretization of the conservation law is then based on 
Stokes formula applied 

div„ g2F= — f^Fx + -^{F^ cos 4,)) . 

cos (p \aX dtp J 

The flux are computed by explicit integration along longitude curves, or latitude curves. 
Length and areas are computed exactly, which allows one to derive a discrete version of 
the geometric compatibility condition. 

5. Perspectives and open problems. In conclusion, based on our theo- 
retical investigations (reported in this brief review), a shock-capturing numerical method 
was designed and implemented in [5j, i.e. a fully intrinsic version of the second-order 
Godunov-type, finite volume scheme. This version, in its principle at first-order accu- 
racy at least, extends to systems like the shallow water equations on the sphere with 
topography. Furthermore, the associated discrete properties (conservation, geometric 
compatibility, asymptotic-preserving) established in |5] extend to systems as well. 

In the setup advocated here for the design numerical schemes, geometric terms in 
hyperbolic conservation laws are taken into account in a direct way by discretizing the 
geometric (weak) formulation of the conservation law, instead of using local coordinates 
and regarding the hyperbolic equation as an equation with variables coefficients and 
source-terms, or instead of viewing the manifold as embedded in a higher dimensional 
Euclidian space. This fully intrinsic approach was validated, both theoretically and 
numerically. It will be interesting to further investigate the properties of these finite 
volume schemes on general curved geometries, especially in the context of complex fiows 
and geometries and to study the large-time asymptotics of solutions. 

In the shallow water equations on the sphere, one unknown — the mass density — 
is a scalar field and can be treated as explained in the present work. However, the 
other unknown — the velocity — is a vector field, and it remains a challenging open 
problem to carry out the general Ricmann problem methodology in this vector- or, more 



HYPERBOLIC CONSERVATION LAWS ON SPACETIMES 



9 



generally, tensor- valued context. This is an important issue, since understanding large- 
scale atmospheric and oceanic motions depends upon designing robust and accurate 
numerical techniques that, we advocate, should directly incorporate curved geometric 
effects. 

The class of geometry-compatible conservation laws on spacetimes forms an analogue 
of the inviscid Burgers equation, which has served as an important simplified model for 
the development of shock-capturing methods. The treatment of boundary conditions on 
non-compact manifolds is an important issue, tackled in |13| , while error estimates a la 
Kuznetsov were derived in l3l I18| . 
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